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ABSTRACT 

We present a comprehensive study of the silicate features at 9.7 and 18 pm of a sample of almost 
800 active galactic nuclei (AGN) with available spectra from the Spitzer InfraRed Spectrograph (IRS). 

We measure the strength of the silicate feature at 9.7 pm, S9.7, before and after subtracting the host 
galaxy emission from the IRS spectra. The numbers of type 1 and 2 AGN with the feature in emission 
increase by 20 and 50%, respectively, once the host galaxy is removed, while 35% of objects with the 
feature originally in absorption exhibit it in even deeper absorption. The peak of S9.7, A pea k, has a 
bimodal distribution when the feature is in emission, with about 65% of the cases showing A pea k >10.2 
pm. Silicates can appear in emission in objects with mid-infrared (MIR) luminosity spanning over six 
orders of magnitude. The derived distributions of the strength of the silicate features at 9.7 and 18 
pm provide a solid test bed for modeling the dust distribution in AGN. Clumpiness is needed in order 
to produce absorption features in unobscured AGN and can also cause the silicates to be in absorption 
at 9.7 pm and in emission at 18 pm in type 1 sources. We find the ‘cosmic’ silicates of Ossenkopf 
et al. to be more consistent with the observations than Draine’s ‘astronomical’ silicates. Finally, we 
discuss the possibility of a foreground absorber to explain the deep silicate absorption features in the 
MIR spectra of some type 2 AGN. 

Subject headings: galaxies: active — infrared: galaxies 


1. INTRODUCTION 


The InfraRed Spectrograph (IRS; [Houck et al. 2004) 


onboard the Spitzer Space telescope provided, over a pe 
riod of more than five years, low- and high-resolution 
mid-infrared (MIR) spectra from many thousands of 
galactic and extragalactic sources, at wavelengths be¬ 
tween 5 and 40 pm. 

In active galactic nuclei (AGN), the MIR emission is 
believed to be UV light reprocessed by the hot dust sur¬ 
rounding the AGN. The dust, often assumed to form a 
toroidal structure in parsec scales around the nucleus in 
a plain extending that of the accretion disk , is consid- 
ered to be distributed either smoothly (e.g. Pier & Kro- 
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Dust is believed to consist of graphite and silicate grams, 
each leaving their unmistakable signature in the Spectral 
Energy Distribution (SED) of AGN, namely the ^1500K 
black-body like rise o f the MIR continuum of typ e (un¬ 
obscured) AGN (e.g. [Hatziminaoglou et al. 2005|), cor¬ 
responding to the sublimation temperature of graphites, 
and an absorption feature centred at 9.7 pm, long known 
to appear in type 2 (obscured) AGN, attributed to sili¬ 
cate grains. 

The uncertain observational evidence for silicate in 


emission in type 1 AGN (Clavel et al. 2000), however, 
posed a problem for the Unified Scheme according to 
which the various types of AGN can be explained by 
alignment effects between the central so urces, the obscur¬ 
ing material (torus) and the observer (Antonucci 1993). 
The problem was finally solved when silicates were un¬ 
ambiguously observed in emission in the MIR spectra of 


many known AGN with IRS 
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Since then, various studies of the behaviour of the sil¬ 
icates have appeared. The v ery first such work s , based 


on a few tens of AGN (e.g. Spoon et al. 2007 Hao et 


M|2007l |Wu et al.|[2009j ) 

feature shows a wide di 


demonstrated that The silicate 
'iversity. On average spectra it 
varies with AGN type, ranging from moderate emission 
in bright quasars to almost no emission or slight absorp¬ 
tion in Seyfert 1 galaxies to stronger absorption in Seyfert 
2 galaxies. Meanwhile, sparse reports on the detection 
of sili cates in emission ( Mason et al.||2009 Nikutta et al. 
2009) showed the diversity of the dust properties, albeit 
some of them rather rare. Several of the above obser¬ 
vations also revealed a second silicate feature at 18 pm, 
as predicted by the models. The relative strength of the 
silicate features at 9.7 and 18 pm has been put forward 
as a possible diagnostic of the torus m orphology (e.g. 
Thom pson et al. 2009j |Feltre et al.||2012|) and chemistry 
(|Sirocky et al.||2008p. 

In this paper, we put together the largest sample of 
active galaxies with available IRS spectra ever composed 
(Sec. [2J) with the aim to complement and extend previ¬ 
ous studies on the MIR characteristics of AGN. To this 
aim, we apply a new spectral decomposition technique 
to separate the nuclear emission from that of the host 
(Sec. [3]). We then proceed with a thorough investigation 
of the"behaviour of the silicate feature at 9.7 pm and 18 
pm in the various AGN types (Sec. [4|. Section [5] dis¬ 
cusses our most important results ana places them into 
a more general context. 
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TABLE 1 

Number of objects per AGN type. 


Type 

N 0 b.i 

Type 

N 0 b.i 

Type 1 AGN 

363 

Syl.2 

24 

Type 2 AGN 

325 

Syl.5 

32 



Syl.8 

18 

SDSS quasars 

141 

Syl.9 

22 


2. THE SAMPLE 


The Cornell Atla S of Spzfcer/Infrared Sp ectrograph 
project (CASSIS^] |Lebouteiller et al. 201 l|) has made 
available the reduced spectra of all the sources observed 
with the low resolution modules of IRS, a total of about 
11000 unique observations. Our master sample is de¬ 
rived from the CASSIS version 6 catalogue, keeping 
each object that fulfills the following requirements: i) 
has an identification in the NASA/IPAC Extragalactic 
Databas^] (NED); ii) has a robust (optical or infrared) 
spectroscopic redshift; iii) the IRS spectrum fully covers 
the range between 6 and 13 micron restframe; iv) the 
median signal-to-noise ratio (SNR) per pixel of the IRS 
spectrum is > 2 . 

To verify requirement i) we rely on the source cross¬ 
identification from CASSIS, which matches the source 
coordinates with the N ED and SIMBAD databases 
(|Lebouteiller et al.||2011|). For sources with no spectro¬ 
scopic redshift in NED, we measure the redshift from the 
IRS spec trum using the template m atching method pre¬ 
sented in Hernan-Caballero (| 2012 |). The typical redshift 
uncertainty with this method is lLz/(l-\-z) ~ 0 . 002 , well 
below the spectral resolution of our resampled spectra 
(AA/A=0.005-0.02). The resampling also increases the 
minimum SNR per resolution element from 2 to > 3 (see 
Sec. |3l ). 


Taking these criteria into account and removing dupli¬ 
cate entries, we end up with a list of 2299 extragalactic 
objects. Out of these, 784 objects have a NED classi¬ 
fication as AGN of types 1, 2 or intermediate, and this 
is the sample we will be working with henceforth. AGN 
of undefined type were left out of the sample. Among 
the type 1 AGN lie 141 quasars from the Sloan Digital 
Sky Survey (SDSS) Data Release 7 Quasar Catalogue 
(Shen et al.| 20 lT|, that for few specific purposes will be 
examined separately. The numbers of the object per sub¬ 
sample are shown in Table [l] The AGN sample is hetero¬ 
geneous but, nevertheless, representative of the infrared 
(IR) AGN population. 


3. SPECTRAL DECOMPOSITION 

The observed MIR spectra of AGN are affected by the 
presence of their host galaxies in two important ways. 
One is the absorption or scattering of AGN emissions by 
material (gas and dust) in the AGN line of sight. This so- 
called foreground absorption modulates the AGN spec¬ 
trum with a multiplicative factor e~ T ^ x \ where r(A) rep¬ 
resents the optical depth at wavelength A. It is not pos¬ 
sible to distinguish foreground absorption from intrinsic 
AGN absorption (that is, the one produced in the AGN 
torus) from MIR data alone, since similar extinction laws 
are considered to apply to dust grains in the torus and 


1 http://cassis.sirtf.com 

2 http://ned.ipac.caltech.edu 


the host. The amount of foreground extinction varies 
from source to source, but at MIR wavelengths it is ex¬ 
pected to be mild in most sources, with the exception of 
some dusty starbursts and edge-on spiral galaxies. The 
other important effect on the AGN spectra is the con¬ 
tamination from host galaxy emission that blends with 
the AGN spectrum. The importance of this background 
emission depends on the relative luminosities of the AGN 
and the host and -crucially- on the spatial resolution of 
the spectroscopic observations. Since the emission from 
the AGN is typically unresolved, an increase in the spa¬ 
tial resolution implies that a larger fraction of the host 
emission can be resolved away. In any case, the back¬ 
ground emission represents an additive modification to 
the AGN spectrum. 

The purpose of our spectral decomposition is to sep¬ 
arate the AGN and host emissions in the integrated 
AGN+host spectra. If successful, this decomposition al¬ 
lows to study the AGN emission as if the host galaxy was 
resolved away. To t his aim we em plo y the de composition 
method presented in Hernan-Caballero et al. (2015). The 
method relies on the large number of high-quality spectra 
in CASSIS to reproduce the spectra of composite sources 
as a linear combination of three CASSIS spectra, each se¬ 
lected from subsamples of sources whose mid-IR emission 
is completely dominated by the AGN, the star-formation, 
or the stellar population. We select these ‘single-spectral- 
component ’ templates as follows: 

For the ‘stellar’ templates, we select 19 local ellipti¬ 
cal and SO galaxies. To ensure they have negligible star 
formation, we require the Polycyclic Aromatic Hydrocar¬ 
bons (PAH) bands to be very weak or absent, with equiv¬ 
alent widths for the 6.2 pm (EW 62 ) and 11.3 pm (EW113) 
PAH features <0.02 pm. We also check that the IRS 
spectra have a blue stellar-like MIR continuum and the 
sources are not classified as AGN in NED. The 54 star¬ 
forming templates (‘PAH’ templates) are IRS spectra of 
normal star-forming and starburst galaxies at redshifts 
up to z=0.14. We make sure these sources do not have 
significant stellar contributions to their MIR spectra by 
requiring both high EW of the PAH features (EW 62 > 
1.0 pm and EW 113 > 1-0 pm) and a very weak contin¬ 
uum at 5 pm. We also verify that they are not classified 
as AGN in NED. Finally, the 147 ‘AGN’ templates are 
IRS spectra of sources classified in the optical as quasars, 
Seyfert galaxies, LINERs, and blazars. We also include 
a variety of optically obscured AGN and radiogalaxies. 
The templates include sources at redshifts from z=0.002 
to z= 1.4 and cover several orders of magnitude in bolo- 
metric luminosity. We ensure that the AGN templates do 
not contain any significant emission from the host galaxy 
by requiring the PAH features to be extremely weak or 
absent (EW 62 < 0.02 pm and EW113 < 0.02 pm). 

Because the AGN and PAH templates are real spec¬ 
tra, each of them already includes some amount of fore¬ 
ground extinction built in. Therefore, we rely on the 
large number of AGN and PAH templates to reproduce 
the diversity of observed spectra that arises from dif¬ 
ferent levels of foreground extinction as well as source 
to source variation in the intrinsic AGN and host spec¬ 
tra. This approach has the advantage of not depend¬ 
ing on assumptions about the -unobservable- intrinsic 
AGN spectrum or the extinction law. Obtaining a good 
fit with this decomposition method requires finding an 
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Fig. 1.— Examples of best-fitting decomposition models for the IRS spectra. The black solid line with grey shading represents the IRS 
spectrum (resampled at AA=0.1 pm resolution) and its l-a uncertainty (photometric errors only). The dotted, dashed, and dot-dashed 
lines represent, respectively, the PAH, AGN, and stellar components of the best-fitting model, shown in yellow. The thin solid line at the 
bottom of each plot represents the residual (spectrum - model). 


AGN template with the appropriate level of foreground 
absorption. This can be problematic for sources with 
very deep absorption features, since few pure-AGN spec¬ 
tra have them. Accordingly, decomposition results for 
sources with deep absorption features have larger resid¬ 
uals and uncertainties. 

We separate the AGN sample into two groups depend¬ 
ing on the spectral coverage: those objects with silicate 
features observed both at 9.7 and at 18 pm in their full 
extent and those for which only the feature at 9.7 pm is 


covered. For sources in the first group, we fit the spectral 
range between 5.2 and 22 pm restframe, while for those 
in the second, we fit only the 5.2 to 15.8 pm range. We 
resample both the spectra and templates to a common 
wavelength grid with an uniform wavelength resolution of 
AA=0.1 pm. This increases the SNR per resolution ele¬ 
ment by ^60% on average, while it still allows to resolve 
important features such as the PAH bands. For every 
galaxy in the sample we try spectral decompositions us¬ 
ing every possible combination of a stellar template, a 
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PAH template, and an AGN template. The best fitting 
model is the one that produces the absolute minimum 
of x 2 . However, to calculate expected values for observ¬ 
ables (e.g. the luminosity of the AGN component or the 
strength of the silicate feature) and their uncertainties, 
we use the full probability distribution functions (PDFs) 
calculated with the ‘max’ method descr ibed i n INoll etl 


ah] ([2009) (for details see |Hernan-Caballero et al.p015, 
She. 2 ). The method also yields, for each object, the 
fractional contribution of each the three components to 
the total luminosity, in the wavelength range of interest. 

Thanks to the use of large sets of real spectra as tem¬ 
plates, our decomposition method manages to reproduce 
the MIR spectrum of composite sources with unprece¬ 
dented accuracy (see Fig. [l]). Typical y 2 values are lower 
than 2 , indicating that residuals in the model fits are 
dominated by noise in both the spectra and templates 
for most sources. 

4. THE SILICATE FEATURES 

The silicate features at 9.7 and 18 pm observed in the 
IR spectra of AGN are believed to arise from the inner, 
hotter parts of the torus or the hot, illuminated side of 
the clumps. Wed efine th e stre ngth of the silicate feature 
following Pier & Krolik ([1992): 


S A =ln 


F(^pea/c) 

Fc(^pea/c) 


(1) 


where F(A pea k) and F c (A pea k) are the flux densities of 
the spectrum and the underlying continuum at the peak 
wavelength of the features, A pea k- A negative (positive) 
value indicates a feature in absorption (emission). 

4.1. The Silicate feature at 9.7 micron 

The top panel of Fig. [2] shows the distribution of the 
strength of the silicate feature at 9.7 pm as measured on 
the original IRS spectra, Sg. 7 tot 5 for type 1 and type 2 
AGN (the 96 AGN of intermediate type are not included 
here and will be discussed separately). 

As already observationally established in the past ten 
years , S 9.7 to t takes a wide range of values (e.g. |Hao et aF 

2007 |Levenso n et al. 2007; [Spoon et ahf 2007 ; Wu et al 


2009, just to name a few). The average spectra of type 

1 AGN exhibit the feature in weak to moderate emission 


(see e.g. Hao et al. 2007 Wu et al. 2009), whil e those 
of type 2 AGN present the feature m absorption (|Sturm 
et al. |2006 ISchweitzer et ah||2008l Mason et al!f [2009 ; 


Hernan-Caballero & Hatziminaoglou 2011). Individually, 

however, type 1 and type 2 AGIN can show silicates in 
absorption and emission, respectively. In our sample of 
698 type 1 and 2 AGN, 35% of the type 1 AGN present 
the feature in absorption and about 15% of the type 2 
AGN show the feature in emission. 

At the same time, when in emission, the peak of the 
feature is often shifted to wavelengths longer than 9.7 
pm in the rest frame, while the shift affects much less 
the featur e when in absorption, as already reported by 
e.gjShi et al.j (|2014|). Fig. |3|shows the distribution of the 
shift, AA pea k = A P e a k — 9.7] pm, as a function of the frac¬ 
tional contribution of the AGN to the luminosity in the 
range between 5 and 15 pm, /agn, for type 1 and type 
2 objects (filled and open symbols, respectively), colour- 
coded by S 9.7 tot- Looking at the sample as a whole, 65% 



S9.7 tot 



Fig. 2.— S 9.7 distribution per for type 1 (blue solid histogram) 
and type 2 (red dashed histogram), before (top panel) and after 
(bottom panel) the subtraction of the host galaxy. 


( 20 %) of the objects with the silicates in emission have 
their A pea k > 10.2 pm (A pea k > 10.6 pm), while the frac¬ 
tion of objects with the same amount of shift among the 
AGN with silicates in absorption is less than 3%. The 
shift to longer wavelenghts is largely associated to a sili¬ 
cate feature in emission, and this in turn only occurs in 
strongly AGN-dominated spectra. 

4.2. Removing the effects of the host 

As S 9.7 tot are measured on the original IRS spectra, 
we expect the derived values to be contaminated by the 
emission of the host galaxy for all objects but those for 
which the AGN completely dominates the MIR emission. 
The top panel of Fig. [4] shows S 9.7 to t as a function of 
/agn- The filled and open symbols correspond to type 1 
and type 2 AGN, respectively. Error bars for S 9.7 tot 
are shown in this plot but they will not be repeated 
in following figures, in order to keep the plots as little 
crowded as possible. What we see here is that as the 
contribution of the host becomes more important (i.e. 
as /agn decreases) the strength of the silicate feature 
decreases, with only few objects exhibiting S 9.7 tot > 0-0 
for /agn <0.7. The dashed line shows the (weak) trend 
for the full sample, with a linear correlation coefficient of 
r=0.45. The trend, however, is driven by type 1 objects 
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Fig. 3.— AA pea k as a function of /agn> for type 1 and type 2 
objects (filled and open symbols, respectively). The symbols are 
colour-coded based on the value of S 9.7 tot- The quantisation of 
A Apeak is an artifact of the algorithm that measures A pea k- 


(filled symbols and corresponding solid line) due to the 
contamination by the emission of the host galaxy. 

In order to see how the silicates behave in the vicin¬ 
ity of the nucleus, we need to remove the contribution 
of the host galaxy, applying the spectral decomposition 
procedure described in Sec. [3j The distribution of the 
strength of the silicate feature on the host-subtracted 
spectrum, S 9.7 agn, is shown in the bottom panel of Fig. 
|2j The behaviour of S 9.7 agn with /agn shown in the 
lower panel of Fig. [4] differs from that of S 9.7 to t in that 
there is now no correlation between the two quantities, 
confirming that the correlation was due to the contami¬ 
nation from the emission of the host, indeed. 

A direct comparison of the two measurements of the 
silicate feature at 9.7 pm is shown in Fig. [5j colour-coded 
by the value of /agn- Objects with MIR emission com¬ 
pletely dominated by the AGN (light-colour symbols) are 
not affected by the subtraction of the host (they lie on 
or very near the 1:1 line). However, as the contribution 
of the host galaxy becomes more important, i.e. /agn 
decreases, (symbols become darker in Fig. [5|, the points 
deviate more and more from the 1:1 line. 

By subtracting the emission of the galaxy, the num¬ 
ber of type 1 AGN with silicates in emission increases by 
20%, reaching 80% of all type 1 AGN, while the num¬ 
ber of type 2 AGN with the feature in emission doubles, 
reaching a total of 25%. At the same time, 35% of both 
type 1 and type 2 AGN with the feature in absorption 
exhibit the feature in even deeper absorption once the 
emission from the host is removed. This happens be¬ 
cause AGN with a silicate feature in deeper absorption 
than that of the surrounding host get their silicate fea¬ 
ture ‘refilled’ in the integrated spectrum. 

In order to check whether S 9.7 agn is af fected by_the 
AGN luminosity, as proposed by e.g. |Maiolino et al. 


(2007), we have to rely on the luminosity at 7 pm, L 7 , 
that spans over six orders of magnitude in our sam¬ 
ple. MIR luminosity in AGN has, in fact, been shown 
to tightly correlate with the X-ray luminosity (see e.g . 
Lutz et aT||2004 Horst et al. 2008 Mateos et al.||2015 ), 



0 0.2 0.4 0.6 0.8 1.0 
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Fig. 4.— S 9.7 as a function of /agn> before (upper panel) and 
after (lower panel) the subtraction of the host galaxy. Filled and 
open circles denote type 1 and type 2 AGN, respectively. The thin 
dashed lines show S 9 . 7 = 0.0 and /agn=0.7, the thick dashed line 
shows the weak (r=0.45) linear correlation of the full sample, and 
the solid and long-dashed lines show the correlations for the type 
1 (r=0.51) and type 2 (r=0.26) objects, respectively. 
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Fig. 5.— S 9.7 tot versus S 9.7 agn, colour-coded based on/AGN- 
Filled and open symbols denote type 1 and type 2 AGN, respec¬ 
tively. 
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S9.7AGN 

Fig. 6.— The relation between S 9.7 agn and L 7 , with points 
colour-coded based on their redshift. The filled and open circles 
correspond to type 1 and type 2 AGN, respectively. 


Fig. 7 . — The distribution of S 9.7 agn for the four intermediate 
types of Seyfert galaxies. 


the ratio of the MIR to the bolometric luminosity de¬ 
pends, however, on the co lumn density alo ng the line 
of sight (Fig. 21 in Hatziminaoglou et al. 2009 ). Fig. 
[ 6 ] shows S9.7 agn as a function of L7 measured on the 
galaxy-subtracted spectrum, with the points coloured 
based on the redshift. There is clearly no dependence 
of S9.7 agn with L 7 and the features can be in emission 
(S9.7 > 0 . 0 ) even in the faintest AGN. Deep silicate fea¬ 
tures (S9.7 <-2) are only found at intermediate luminosi¬ 
ties (or redshifts). A deeply obscured but low-luminosity 
AGN would be overwhelmed by the emission of its host, 
that would ‘refill’ the silicate feature. The lack of deep 
silicates in objects with very high IR luminosities and/or 
high redshifts, on the other hand, suggests a selection 
effect, as they might be too faint in the optical for a 
reliable identification. 


4 . 2 . 1 . AGN of intermediate types 

Among the 784 AGN of our sample, 96 have a NED 
classification of intermediate type Seyfert galaxies. The 
number of objects per type (comparable in all four sub¬ 
samples) is shown in Table [lj Figure [7] shows the distri¬ 
bution of S 9.7 agn for the four sub-samples. 

The mean average S 9.7 agn decreases from Syl.2 to¬ 
wards later types, with Sy 1.2, 1.5 and 1.8 showing on 
average the silicate feature in very weak or no emission, 
while Syl.9 has a negative average value. Additionally, 
as we move to later intermediate types, the dominance 
of the AGN component to the MIR emission decreases: 
while 96% of Syl.2 objects have /agn > 0.7, the frac¬ 
tion dr ops to 55% for Syl.9s. This is in agreement with 


Deo et al. (2007) that found the MIR spectra of Syl .8 
and by 1.9 to be dominated by starburst features (PAH). 
The mean values of S 9.7 agn and the fraction of objects 
with /agn > 0.7 for each intermediate type are shown 
in Table HI Note that the shift towards lower values of 
S 9.7 agn when moving towards later Seyfert types persist 
even when only objects with /agn >0.7 are considered, 
as seen in the right-most column of table [ 2 ] Finally, 
the Ap ea k of intermediate Seyfert types shows the same 
behaviour as for the rest of the AGN sample, i.e. as 


TABLE 2 

Mean values and standard deviations of S 9.7 for each 

INTERMEDIATE AGN TYPE, FRACTION OF OBJECTS WITH 
/AGN >0.7 AND S 9.7 FOR THAT FRACTION. 


Type 

(S 9.7 tot) 

/agn > 0-7 

/c/agn> u - 7 \ 

W 9.7 tot / 

Syl .2 

0.175T0.13 

96% 

0.160T0.14 

Syl.5 

0.043T0.46 

90% 

0.087T0.22 

Syl .8 

0 . 022 T 0.22 

61% 

0.059T0.il 

Syl.9 

-0.305T0.49 

55% 

-0.227T0.35 


described in Sec. 14.11 


4 . 2 . 2 . SDSS quasars 

Out of the 784 AGN, 141 are spectroscopically con¬ 
firmed SDSS quasars, for which estimates of the mass of 
the central black hole, Mbh, derived from emission line 
measurements, as well as bolometric lumino sities, Lboi, 


deriv e d from fitting techni ques are available (Shen et al. 
2011). Maiolino et al. (2007) reported an increase of 89.7 
with increasing IVIbh, from low-luminosity, low-redshift 
type 1 AGN to high-lu minosity, high-redshift quasars. 
Thompson et al. (2009), on the other hand, found no 
trend of with luminosity. Mbh for the 141 SDSS quasars 
in question spans the range [I Q 7,3 M^ — 10 1Q M^ j, i.e. 


Maiolino et al. 


(2007) sam 


almost identical to that of the ! 
pie, but we do not find any correlation of"S 9 . 7 agn with 
Mbh, as shown in Fig. [ 8 ] Note, however, that the bolo¬ 
metric luminosities of the SDSS quasars of our sample 
are all above 8.2 xlO 10 L e (10 44 , 5 erg/sec), i.e. the two 
samples are not directly comparable. In the quasar sub¬ 
sample, S 9.7 agn and Lboi are completely unc orrelated 
(see colour-co ding in Fig. [ 8 |), in agreement with Thomp¬ 
son et al. (2009). 


4 . 3 . The Silicate feature at 18 micron 


The silicate feature at 18 pm, though predicted by 
models and observed in many of the low-to-intermed iate 


redshift AGN (z typically < 0.5) by now (e.g. 

[Hao et 

al. 2005 

Sirocky et al. 2008 Thompson et al. 20( 

39|, has 


received less attention than its lower wavelength coun- 
terpart, as its measurement presents a greater challenge. 
For one, it often overlaps with the MIR bump of the 
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M B h [M,] 

Fig. 8 .— S 9.7 agn as a function of Mbh for the 141 SDSS quasars 
of the sample. The points are colour-coded by Lbol- 



Fig. 9.— Sis agn distribution per AGN type (type 1 in blue 
solid histogram, type 2 in red dashed histogram). 


AGN SED (e.g. Prieto et al. 


2010 ). Also, the steep mid- 


to-far IR emission from the host implies it contributes a 
higher fraction to the continuum emission at these longer 
wavelengths, making the measurement of the feature a 
tedious job. Following the procedure described in Sec. 
[3j but extending the wavelength range to 22 pm rest- 
frame, we perform spectral decompositions and measure 
the strength of the 18 pm silicate feature in the AGN 
component, Sis agn> for the 631 AGN with adequate 
wavelength coverage. Fig. [9] shows the distribution of 
Sis agn for type 1 (blue solid histograms) and type 2 
(red dashed histograms). 

While at least as prominent as its counterpart at 9.7 
pm in emission, when in absorption the feature only 
reaches moderate depths. Furthermore, more than 50% 
of type 2 AGN exhibit the feature in emission while only 
10% of type 1 AGN have it in absorption. Comparing 
these numbers with those for S 9.7 agn it becomes obvi¬ 
ous that the feature can be in absorption at 9.7 pm while 


still in emission at 18 pm. We will get back to this point 
in Sec. [3 


5. DISCUSSION AND CONCLUSIONS 


Using the largest sample of AGN with MIR spec¬ 
troscopy ever assembled we quantify, for the first time, 
the effects of the emission of the host galaxy on the be¬ 
haviour of the silicate features, a tracer of the properties 
of the hot dust in the torus. We rely on the classifica¬ 
tion provided by NED in order to call an object “AGN” 
as well as for the classification of AGN in types 1 and 
2 . The sample includes a variety of AGN, from objects 
where the AGN completely dominates the MIR emission 
(/agn ^ 1) to AGN whose MIR emission is almost en¬ 
tirely dominated by star formation (/agn << !)• We 
find the fraction of objects with a strong contamination 
from the galaxy (/agn < 0-7) to be much higher among 
type 2 AGN (43 pre cent) than among type 1 AGN (12 
pre cent). 

The emission of the host affects the behaviour of the sil¬ 
icate features. Broadly speaking, the strength of the sil¬ 
icate feature at 9.7 pm is a measure of the optical depth, 
T 9 . 7 , along the line of sight, and goes from em ission to 


absorption as T 9.7 increases (see e.g. Fritz et al. 2006 
Fig. 9). In this simple picture, type 1 and type 2 AGbl 
should show the feature in emission and absorption, re¬ 
spectively. Our study, however, shows that type 1 (type 
2 ) AGN with the feature in absorption (emission) are 
very common, even after subtracting the contribution of 
the host galaxy. The numbers of type 1 and 2 AGN 
with the feature in emission increase by 20 and 50%, re¬ 
spectively, once the host galaxy is removed, while 35% 
of objects with the feature originally in absorption ex¬ 
hibit it in even deeper absorption after subtraction of the 
host. This means that the combined spectrum exhibits 
an S 9.7 tot intermediate between S 9.7 agn and that of the 
host. The host galaxy nearly always shows mild silicate 
absorption, as the power sources (stars) are well mixed 
with the absorbing dust, and therefore the strength of 
the silicate feature does not correlate with the optical 
depth of the gas and dust along the line of sight. Conse¬ 
quently, contamination from the host increases the depth 
of the absorption if the feature intrinsic to the AGN is 
in emission or mild absorption, but decreases the depth 
(i.e. fills the gap) if the AGN has a deeper feature than 
the host. 


S9.7 tot is scarcely ever shown in emission when the 
MIR emission is strongly contaminated by the host 
galaxy (/agn < 0-7), with those objects exhibiting on 
average the feature in deeper absorption than their AGN- 
dominated counterparts. S 9 . 7 AGN, on the other hand, 
shows no dependency on /agn and appears from moder¬ 
ate emission to quite deep absorption with all possibili¬ 
ties in between. The lack of correlation between S9.7 agn 
and /agn 5 two quantities that are physically unrelated, 
indicates that the decomposition mechanism successfully 
removes the most important part of the host galaxy emis¬ 
sion, allowing for an almost unbiased silicate measure¬ 
ment of the silicate feature in the vicinity of the nucleus. 

We have also addressed the issue of the observed shift 
of the silicate feature at 9.7 pm to longer wavelengths 
with respect to the nominal peak at 9.7 pm. Our decom¬ 
position method does not allow reliable estimate of A pea k 
on the galaxy-subtracted spectra. We therefore rely on 























Hatziminaoglou et al. 


the measurement carried out on the original IRS spec¬ 
tra. We find the largest shifts (A pea k > 10.2 pm) to ap¬ 
pear only in objects with an important AGN component 
(/agn >0.7) with the feature in emission, regardless of 
their type. When the feature is in absorption, and again 
irrespective of the type, it appears at or near its nomi¬ 
nal wavelength (A pea k < 10.2 pm). Since its discovery, 
various scenarii have been proposed in order to explain 
this shift, such as the presen ce of porous dust (Li, Shi 
& Li||2008| iS mith et al.|| 2010 |), the presenc e of different 
dust species (Markwick-Kemper et al. 2007 ), or radiativ e 
transfer effects (IJNikutta et al. 2009|. Shi et al. (2014), 
however, posit that this is not a radiative transfer effect 
but rather the effect of direct exposure of the silicates to 
the nuclear radiation, that modifies the size or the chem¬ 
ical composi tion of the grains, more along the lines of 
Smith et al.l (120101). 


5.1. Dust distribution inside and outside the torus 

S 9.7 agn is seen in only moderate emission or slight 
absorption in most type 1 AGN, a behaviour that tra¬ 
ditionally favours a clumpy morphology, both because 
smooth models often predict silicates in stronger emis¬ 
sion than observed and because, with the exception of 
a couple of smooth model parameter combinations, only 
clumpiness can explain the feature in absorption in un¬ 
obscured AGN. However, and even though not produc¬ 
ing the feature in absorption, a continuous dust distribu¬ 
tion can also give rise to silicates in only weak emission 
for a large variety of parameters, as shown in Fi g . 4 of 
Feltre et al. ( 2012 ). Furth ermore, Sirocky et al. (2008) 
showed that the use of the Ossenkopf et al.| Q1992 ) dust 
absorption and scattering coefficients result in consider¬ 
ably less prominent silic ate emi ssion features compared 
to other dust models like Draine (2003). We therefore in¬ 
terpret the behaviour o f $ 9.7 ag n alone as favo uring the 
Ossenkopf et al. (1992) over the Draine (2003) silicates, 
but this property alone cannot provide much insight into 
the morphology of the dust. 

The transition of the mean values of S 9.7 from weak 
emission to absorption from Syl .2 to Syl.9 can be 
explained by either dust morphologies: in a smooth 
medium it can be attributed to an increase of the in¬ 
clination (as measured from the poles) and hence the 
intervening material, with the silicate feature at 9.7 pm 
that arises from the inner, hotter parts of the torus, being 
increasingly blocked by the bulk of the dust as the view¬ 
ing angle increases. In a clumpy medium, on the other 
hand, this could simply be attributed to different levels 
of obscuration along the line of sight, independently on 
the orientation. 

The combined strength of the silicate features at 9.7 
and 18 pm is sensitive to the chemistry and mor phology 
of the dust surrounding the AGN, i.e. the torus dSir ocky 


et al.|2008||Thompson et al.|2009 |Feltre et al.| 2012 |). To 
test this, we created th ree grids of models, t wo c lumpy 
and a smoot h, following Nenkova et al. (2008) and Feltre 
et al. (2012), respectively. All three sets o f models share 
the same the primary source, described in Nenkova et al. 
(2008). One of the clumpy grids was created using the sil¬ 


i cates absorption and scattering coefficients from Draine 
(2003) while the other c lumpy as well as the smo oth grids 
were created using the Ossenkopf et al. (1992) silicates. 
The models have been created to have matched param- 



J> -0.4 


- 0.6 

- 0.8 

- 1.0 



Fig. 10.— Si 8 agn as a function of S 9.7 agn> overplotted on 
model predictions. The top (bottom) panel shows type 1 (type 
2 ) objects (in blue). Smooth models are shown in green, clumpy 
models with Ossenkopf et al. (1992) and Draine (2003) silicates 
are shown in grey and pink, respectively. 


eters, as defined in Feltre et al. ( 2012 ), i.e. each of the 


smooth models in the grid has an~equivalent model (in 
terms of geometrical properties) in the clumpy grid. The 
parameter space explored by the model grids is briefly 
described in Appendix [Aj 

Fig. [lO] shows the distributions of Sig agn and 
S 9.7 agn 5 compared to model predictions. The top (bot¬ 
tom) panel shows type 1 (type 2 ) AGN (in blue). Smooth 


models are shown in 1 

?reen, clumpy models using the 

Ossenkopf et al. ( 

1992 

| and Draine (|2003) silicates are 

shown in grey anc 

pink, respectively. T he spread of the 


observed data points indicates the variety of torus ge¬ 
ometries in nature in terms of size, shape and optical 
depths. 

As previously noted b y ot he r aut hors, the so called 
‘astronomical’ silicates (Draine 2003) produce stronger 
emission features at 9.7 pm and a wider range of 
Sig AGN/S 9.7 agn than wha t is observed, compared to 
the Ossenkopf et al. (1992) silicates. The bulk of type 
1 objects lie in the region of overlap between smooth 
and clumpy models. Overall, models reproduce in a bet¬ 
ter fashion the silicates in type 2 AGN. However, both 
smooth and clumpy models, albeit covering a large pa- 
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rameter space, fail to reproduce the deep absorption fea¬ 
tures seen in type 1 AGN (lower left quadrant in top 


panel of Fig. \w\j. At the same time, and as already 
suggested by Imanishi et al.| (|2007|), a continuous dust 
distribution is the only morphology that can reproduce 
the deepest absorption features seen in the sample for 
typ e 2 views (lower left quadrant of bottom panel in Fig. 


10 ), at least without resortin g to additional obscur ation 
by the host galaxy. In fact, |Levenson et al.| ([2007) sug¬ 
gested that a deep silicate absorption feature at 9.7 pm 
requires the primary source to be embedded in a contin¬ 
uous, optically and geometrically thick dusty medium, 
while a clumpy medium, with clouds illuminated from 
the outside, will result in a much shallower feature, as 
the emission from these clouds will fill the absorption 
trough. 

We visually inspected all AGN with S 9.7 agn < —1-5 
for which we could find high enough resolution images 
and established that they are predominantly hosted in 
galaxies with high inclinations, galaxies with very vis¬ 
ible dust lanes crossing the centre (e.g. NGC5793 or 

NGC7172), or reside in interactive systems at various 
stages of their inter action (e.g. Mkr 273 , Mrk 331, 
NGC2623). In fact, Goulding et al. ( 2012 |) reach this 
same conclusion in their study of 20 nearby (z < 0.05) 
bona-fide Compton-thick AGN. This implies that the 
source of the deepest silicate absorption features is dust 
in the host galax y rather than dust in the t orus (see also 


Deo et al. 


2007). In favour of this view, Lagos et al 


(| 201 l|) find that type 1 AGN have a strong preference in 
residing in face-on galaxies, while the type 2 AGN reside 
in hosts of any orientation. The deep absorption fea¬ 
tures are, therefore, of no relevance to the modeling of 
the torus and they should not be used to favour smooth 
dust distributions over clumpy ones. We note that this 
is by no means an argument against the AGN unifica¬ 
tion scheme, but does suggest that the obscuration of 
the nucleus, i.e. a type 2 view, may, in some cases, be 
decoupled from the orientation of the torus. 

Finally, what is not reproduce by none of the torus 
models are many of the obscured and unobscured AGN 
with the silicates absorption at 9.7 pm but in emission at 
18 pm. Smooth models do not predict at all such a be¬ 
haviour in type 1 views, while both morphologies produce 
Si8 agn of about half the strength than that measu red 
on the spectra when S 9.7 agn is in absorption. F elt re et| 


al. ( 2012 ) showed that the adopted primary source can 
also affect the properties of the silicate feature. As a last 
resort, and since we have no control of the primary source 
of the Nenkova models, we produced an addition al grid of 
smooth models with the |Ossenkopf et al.|([1992|) s ilicates 

3 et al.l ([2 


but using the primary source from Felt re et al. (|2012), 
that allows for more intrinsic AGN emission at wave¬ 
lengths beyond 1 pm c om pared to the primary source 
used by the Nenkova et al. (2008) models, which we have 
adopted for the other grids. These models still fail to re¬ 
produce the objects in the upper left quadrant of the fig¬ 


ure for type 1 views, they do however reproduce a much 
larger fraction of the objects in that same quadrant for 
type 2 views than any of the other model grids presented 
before, as shown in Fig. [Tl] 

To summarise, and even though the emerging picture 



Fig. 11.— Si8 agn as a function of S 9.7 agn for type 2 objects 
(in blue), overplotted on smooth model predictions (in yellow), 
using the Feltre et al. (2012) primary source. 


satisfactory fashion the absorption and scattering prop¬ 
erties of the silicates in the obscuring torus. Clumpi¬ 
ness is needed in order to produce absorption features in 
unobscured AGN, if no foreground absorber is invoked. 
Clumpiness can also cause the silicates to be in absorp¬ 
tion at 9.7 pm and in emission at 18 pm in type 1 sources, 
but a primary source with more intrinsic AGN emission 
at A > 1.0 pm might be necessary to create stronger 
Sis agn in emission, when S 9.7 agn < 0 . 0 . 
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APPENDIX 

SMOOTH AND CLUMPY MODEL PARAMETERS 

The two original smooth and clumpy model grids cover large but not entirely overlapping parameter space. For the 
purposes of this work, we matched the two sets of parameters to produce models that can be directly comparable. 
Both models share the inner-to-outer radius ratio, Y . The radial distribution of the dust or the clumps, g, can also 
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TABLE 3 

Matched SMOOTH AND CLUMPY MODEL PARAMETERS. 



Smooth 

Clumpy 

Y 

10, 30, 60, 100, 150 

10, 30, 60, 100, 150 

q 

0 , 1 

0 , 1 

7 

1.65, 2.25, 3.24, 5.06, 9.0 


(7 


70°,60°,50°,40°, 30° 

T9.7 

0.21, 0.42, 0.63, 0.84, 1.05, 
1.26, 1.47, 1.68, 1.89, 2.1, 
2.52, 2.94, 3.36, 3.78, 4.2, 
5.04, 5.88, 6.3, 6.72, 7.56, 
8.4, 10.08, 11.34, 11.76, 12.6 


N 0 


1-5 

TV 


5, 10, 20, 30, 40, 

60, 80, 100, 150 


be considered as equivalent between the two models. For smooth models, the dust distribution in polar coordinates is 
given by: 

p(r,9) = p 0 • • e- 7X|cos(e)l (Al) 

while clumpy models assume a Gaussian angular distribution of clouds of width a given by: 

N T {p) = N 0 e^^\ (A 2 ) 

where Nq is the average number of clouds along a given radial direction on the equatorial plain and /3 (=90-0) is the 
angle with respect to the axis of the torus. Considering the two extreme viewing angles for both models, the torus 
opening angle is of no relevance and 7 and a can be taken such to match each other. While smooth models use the 
equatorial optical depth at 9.7 pm, 79.7, as an input parameter, clumpy models use the optical depth of a single could 
calcul ated in the W -band at 0.55 um, ry. 79.7 for clumpy models can be derived via the relation 79.7 = 0.042 xryxiVo, 
(see Nenko vaet al.|2 008). Table [3] summarises the values of each of the parameters for the matched model grids. 


REFERENCES 


Antonucci R. 1993, ARA&A, 31, 473 

Buchanan C. L., Gallimore J. F., O’Dea C. P. et al. 2006, AJ, 

132, 401 

Clavel J., Schulz B., Altieri B.et al. 2000, A&A, 357, 839 
Deo R. P., Crenshaw D. M., Kraemer S. B. et al. 2007, ApJ, 671, 
124 

Draine B. T. 2003, ApJ, 598, 1017 

Feltre A., Hatziminaoglou E., Fritz J., Franceschini A. 2012, 
MNRAS, 426, 120 

Fritz J., Franceschini A., Hatziminaoglou E. 2006, MNRAS, 366, 
767 

Goulding A. D., Alexander D. M., Bauer F. E. et al. 2012, ApJ, 
755, 5 

Granato G.L., Danese L. 1994, MNRAS, 268, 235 
Hao L., Weedman D. W. Spoon H. W. W. et al. 2007, ApJL, 655, 
L77 

Hao L., Spoon H. W. W., Sloan G.C. et al. 2005, ApJL, 625, L75 
Hatziminaoglou E., Fritz J., Jarrett T. 2009, MNRAS, 399, 1206 
Hatziminaoglou E., Perez-Fournon I., Polletta M. et al. 2005, AJ, 
129, 1198 

Hernan-Caballero A., Alonso-Herrero A., Hatziminaoglou E. et al 
2015, ApJ, in press 

Hernan-Caballero A. 2012, MNRAS, 427, 816 
Hernan-Caballero A. &; Hatziminaoglou E. 2011, MNRAS, 414, 
500 

Hernan-Caballero A., Perez-Fournon I., Hatziminaoglou E. et al. 
2009, MNRAS, 395, 1695 

Hiner K. D., Canalizo G., Lacy M. et al. 2009, ApJ, 706, 508 
Honig S.F., Beckert T., Ohnaka K. et al. 2006, A&A, 452, 459 
Horst H., Gandhi P., Smette A., Duschl W. J. 2008, A&A, 479, 
389 

Houck, J. R., Roellig, T. L., van Cleve, J., et al. 2004, ApJS, 154, 
18 

Imanishi M., Dudley C. C., Maiolino R. et al. 2007, ApJS, 171, 72 
Lagos C., Padilla N. D., Strauss M. A., Cora S. A., Hao L. 2011, 
MNRAS, 414, 2148 

Lebouteiller V., Barry D.J., Spoon H.W.W. et al. 2011, ApJS, 
196, 8 


Levenson N. A., Sirocky M. M., Hao L. et al. 2007, ApJL, 654, 
L45 

Li M. P., Shi Q. J., Li A. 2008, MNRAS, 391, L49 

Lutz D., Maiolino R., Spoon H. W. W., Moorwood A. F. M. 

2004, A&A, 418, 465 

Maiolino R., Shemmer O., Imanishi M. et al. 2007, A&A, 468, 979 
Markwick-Kemper F., Gallager S. C., Hines D. C., Bouwman J. 
2007, ApJL, 668, L107 

Mason R. E., Levenson N. A., Shi Y. et al. 2009, ApJL, 693, L136 
Mateos S., Carrera F. J., Alonso-Herrero A.et al. 2015, MNRAS, 
in press 

Nenkova M., Sirocky M. M., Ivezic Z., & Elitzur M. 2008, ApJ, 
685, 147 

Nikutta R., Elitzur M., & Lacy M. 2009, ApJ, 707, 1550 
Noll S., Burgarella D., Giovannoli E. et al. 2009, A&A, 507, 1793 
Ossenkopf V., Henning T., & Mathis J. S. 1992, A&A, 261, 567 
Pier E.A. & Krolik J.H. 1992, ApJ, 401, 99 
Prieto M. A., Reunanen J., Tristram K. R. W. et al. 2010, 
MNRAS, 402, 724 

Schweitzer M., Groves B., Netzer H. et al. 2008, ApJ, 679, 101 
Shen Y., Richards G. T., Strauss M. A. et al. 2011, ApJS, 194, 45 
Shi Y., Rieke G. H., Ogle P. M., Su K. Y. L., Balog Z. 2014, 
ApJS, 214, 23 

Shi Y., Rieke G. H., Hines D. C. et al. 2006, ApJ, 653, 127 
Siebenmorgen R., Haas M., Kriigel E., Schulz B. 2005, A&A, 436, 
L5 

Sirocky M. M., Levenson N. A., Elitzur M., Spoon H. W. W., & 
Armus L. 2008, ApJ, 678, 729 
Smith H. A., Li A., Li M. P. et al. 2010, ApJ, 716, 490 
Smith J. D. T., Draine B. T., Dale D. A. et al. 2007, ApJ, 656, 
770 

Spoon H.W.W., Marshall J.A., Houck J.R. et al. 2007, ApJ, 654, 
L49 

Sturm E., Hasinger G., Lehmann I. et al. 2006, ApJ, 642, 81 
Sturm E., Schweitzer M., Lutz D. et al. 2005, ApJL, 629, L21 
Thompson G. D., Levenson N. A., Uddin S. A., & Sirocky M. M. 
2009, ApJ, 697, 182 

Wu Y., Helou G., Armus L. et al. 2010, ApJ, 723, 895 










Silicate features in active galaxies 


11 


Wu Y., Charmandaris V., Huang, J., Spinoglio, L., Tommasin, S. 

2009, ApJ, 701, 658 



